Calculation of the band gap energy of ionic crystals 



00 
On 
OS 

D 

c/3 



A. Aguado 1 , A. Ayuela 2 , J. M. Lopez 1 , J. A. Alonso 1 , J. F. Rivas-Silva 3 and M. Berrondo 4 
1. Departamento de Fisica Teorica, Facultad de Ciencias, Universidad de Valladolid, 47011 Valladolid, Spain. 
2. Institut fur Theoretische Physik, Technische Universitdt Dresden, 01062 Dresden, Germany. 
3. Instituto de Fisica, Universidad Autonoma de Puebla. 72750 Puebla, Mexico. 
4- Brigham Young University, Provo, UT84602, USA. 
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using the Perturbed-Ion model supplemented with some assumptions for the treatment of excited 
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I. INTRODUCTION 

The increasing interest which is presently observed in 
luminescent materials is due to their numerous technolog- 
ical applications, namely, luminescent lighting and prepa- 
ration of lamp phosphors, nuclear spectroscopy, laser sci- 
ence, or the construction of two-dimensional detectors for 
usejn medical screens and crystallography, to mention a 
fewtl. The calculation of the properties of those materi- 
als requires a model for the material and a computational 
model. The first one concerns the modelling of the crys- 
tal of interest. The computational model concerns the 
level of theory used in the calculations. The calculation 
of the electronic properties of doped ionic crystals (or 
crystals containing vacancies) at an ab initio level is still 
a challenge for the computational methods normally used 
for crystalline solids or molecules. The reason in the first 
case is the breaking of the translational symmetry by the 
presence of the defect. This enters in conflict with Bloch's 
theorem, which is the basis of the solid state calculations. 
The use of molecular methods in real space requires a 
large number of atoms to be treated in a selfconsistent 
wayn. Those molecular methods replace the whole crys- 
tal by a cluster, but the role of the surface should not be 
overlooked even for large clusters. An intermediate ap- 
proach is provided by embedded cluster models, in which 
the impurity is surrounded by a small fraction of the crys- 
talline environment, and the rest of the crystal is simply 
described by point charges. This approach has been al- 
ready applied to the problem of luminescent impuritiesH. 
Evidently, the description by point charges only gives an 
approximate representation of the bulk. In summary, a 
large number of atoms is needed for simulating an ionic 
crystal containing an impurity, and even in such a case 
the inherent problems of the cluster surface and the de- 
scription of the infinite bulk should be faced. Thus, a 
model that gives a better description of the impurity- 
lattice and cluster-lattice interactions is desirable. The 
ideas of the Theory of Electronic Segaxability (TES) de- 
veloped by Huzinaga and coworkersaTl should be useful 
in this context. The TES supplies a natural framework 



to develop accurate schemes for dealing with thejpluster- 
lattice interaction. The perturbed ion (PI) modelQ, based 
on the TES, was developed for the study of ionic crystals 
or, more in general, for crystals formed by atoms with 
closed shells, and has ijpgen later extended to free clusters 
of the same materialsBfi 

We plan to apply the PI model to study luminescent 
centers in ionic crystals. As a first step, we use the PI 
model in this paper to calculate the band gap of pure 
alkali halide, alkaline-earth oxide, AI2O3 and Si02 crys- 
tals. All these materials are important in the field of 
luminescence and also in other fields ranging from catal- 
ysis to magneto-optical devices. The calculation of the 
energy band gap of these materials is important for sev- 
eral reasons: first of all, the impurity levels of the doped 
crystals are located in the band gap of the pure crystal. 
Furthermore, the gap is a very important quantity in the 
first step of the scintillation process, namely, absorption 
of radiation leading to formation of electron- hole pairs. 
This step influences the global efficiency of the scintil- 
lator. Last but not the least, it is also important from 
the theoretical point of view, as it gives information on 
the quality of the model. A large body of work exists on 
the calculation of the band gaps in these materials us- 
ing the traditional methods of band theorycTcj, and we 
compare our results to a representative set of band cal- 
culations. Although we obtain a better global agreement 
with experimental gaps than many ab initio band struc- 
ture calculations, our intention here is not to compete 
with those well established solid state methods. Instead, 
we only intend to show that the PI model affords an accu- 
rate description of the band gap within the framework of 
a cluster-like approach, which is considered a convenient 
approach for the study of doped crystals. 

The structure of the paper is as follows: In section II 
we present the PI model for ionic crystals, showing how 
several problems concerning the cluster approximation 
are circunvented. Since the PI model is originally formu- 
lated in a Hartree-Fock(HF) framework, we also discuss 
the introduction of Coulomb correlation. Section III de- 
scribes the calculation of the gap. Results are presented 
for the above mentioned crystals and the trends obtained 



I 



are discussed in comparison with experiment and other 
calculations. Finally, section IV summarizes our conclu- 
sions. 



II. THE PERTURBED ION MODEL 
A. Theory 

The PI model has been developed for systems (pure 
crystals, crystals with defects, or finite clusters) formed 
by weakly correlated one-center electronic groups, the 
prototypical systems being ionic crystals, like the alkali 
halides, formed by closed-shell ions. Therefore, accord- 
ing to the theory of electronic separabilityBcl, the wave- 
function of the system can be expressed as an antisym- 
metrized product of the local wavefunctions describing 
each group. If these loftalpwavefunctions satisfy strong- 
orthogonality conditionsc-IE3, the total energy is the sum 
of intragroup, or net energies, and intergroup, or inter- 
action energies: 

E system = E E% et + 2 E E E *nt ^ 
R R S(=£R) 

where the R and S sums run through local groups. All 
the contributions to E S y S tem due to a group A can be 
collected in an effective energy defined as 

E tff ~ E net + E ^tnt = ^net + ^tnt > (2) 
S&A) 

but it should be noted that E sy stem is not a simple sum 
of the effective energies since the interaction energies are 
then counted twice. So the additive energy of a group A 
is defined as 

Eadd — Enet + 2 E = E ™t + T^tnt (3) 

S&A) 

and then the total energy of a system formed by a groups 
of type A, b groups of type B, etc, can be written 

E system {A a B b ...) = aE add + bE^ dd + ... (4) 

In practice, for ionic crystals the local groups will be 
identified with the closed-shell ions. 

Group wavefunctions are obtained by minimizing their 
effective energies if strong orthogonality conditions are 
satisfied among the group wavefunctions. This restricted 
variational procedure can be successive and iteratively 
applied to all the different groups, in order to determine 
fully consistent group wavefunctions and the best sys- 
tem wavefunction compatible with the initial assumption 
of separability. When the wavefunction of a particular 
group is being determined, the wavefunctions of all the 
other groups are considered frozen. 

The effective energies can be expressed as expectation 
values of appropriate effective Hamiltonians. For group 
A 

Et ff =< i> A | H? ff | Vu > - J2 Z A v£(Ras). (5) 

S&A) 



The group- wavefunction if; a is a Slater determinant. The 
first term in eq.(5), namely < if) a I ^eff I V'A >, collects 
the energy associated to the electronic cloud of group A, 
and the second term gives the interaction between the 
nucleus A, of charge Z A , and the Hartree potential of 
groups A). Ras —\ Ra — Rs | indicates the distance 
between the two nuclei. The effective Hamiltonian H^f 
can be written as a sum of several terms 

N A N a 
i=l t=l 

N A 

E tf + E E + **(<)]. (6) 

l<j<i<N A »=1 S&A) 

with TiA =| ?i — Ra I and ry =| fi — fj \. The i and 
j-summations run over the Na electrons of group A. The 
first three terms in the Hamiltonian represent the elec- 
tronic kinetic energy, the interaction between the elec- 
trons and the nucleus and the interelectronic repulsion 
(those three terms refer exclusively to group A). On the 
other hand, the terms containing the S-summation ac- 
count for the interaction with the frozen groups S. This 
interaction separates itself in two parts. First, VffAi) 
represents the effective potential energy of an electron 
(of group A) in the mean field of the group S: 

V e s ff (i) = -Z s rr s } + vS(i) + V|(i) = V*(i) + y£(i), 

(7) 

where the three terms in Vj^y are the electron-nucleus, 
the classical electron-electron and the exchange parts of 
the potential energy (the sum of the first two terms is the 
Hartree potential V^). Second, the strong orthogonality 
between the orbitals of the active group A and those of 
the other groups S is enforced in H e / / by means of the 
projection operator P (i). For systems-formed by closed 
shell ions, this operator takes the formE£l 

p5 (*) = Ei^>(- 2e 9 5 )<^i< ( g ) 

where g runs over all occupied one-electron orbitals <p g 
(with orbital energies ) of group S. 

We now deal with the explicit form of V^^El. For the 
closed-shell ions considered here the Hartree part of this 
potential is given by 

V§ (fi ) = - Z s r-l +Jp s (r 2 )r 1 - 2 1 dr 2 (9) 

where p s (r) is the electron density of the ion. It is phys- 
ically and computationally convenient to separate this 
electrostatic potential into classical and nonclassical (nc) 
terms such that: 

V$ (ns) = -q s r^ + Vl(r is ) (10) 

where q s is the net charge of the ion, and V£ c represents 
the deviation of V§ from a point charge potential due 
to the finite extension of the electronic density of the 
ion. This operator can be effectively computed using 
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the Sflverstone expansion of a function in a displaced 
centeiOH 

The exchange operator VJL can be written as the non- 
diagonal spectral resolutionEj 

I 

v|(i) = - ^2 J2 Hi alm > s > A ( l > ab ' < blm > 5 1 • 

l m— — l a, 6 

(ii) 

Here | aim, S > are products of spherical harmonics Yj m 
and primitive radial functions for the S ion, a and b run 
over the Slater- Type Orbitals (STO) of I symmetry, and 
A(l, ab, S) are the elements of the matrix 

A = S- 1 K.S- 1 (12) 

where S and /C are the overlap and the exchange matrices 
for the S ion in the {| aim, S >} basis. Full details-,can 
be found in the original paper by Luaha and PueyoH. 

As indicated above, when the effective Hamiltonian 
of an active group A is being diagonalized, all the 
other groups S are considered frozen. Selfconsistency is 
achieved following an iterative scheme. At each itera- 
tion all the inequivalent ions are successively treated as 
active groups, and the orbital wavefunctions and eigen- 
values of the frozen groups are taken from the previous 
iteration. The iterations are continued until convergence 
is achieved. 

Some features of the PI model are worth mentioning 
here. In the cluster approaches the solid is partitioned 
into a "cluster region" and "the rest of the crystal" . The 
physical and mathematical description of both entities 
should be accurate enough, as well as their interrela- 
tions. Frecuently, the cluster has been solved using pre- 
cise quantum-mechanical methods while therenvironment 
has been simulated by using point chargescj, but more 
accurate descriptions of the environment surrounding the 
active cluster are necessary in-erder to account properly 
for cluster-lattice interactional. Besides, the election of 
the cluster size is also a delicate problem, because surface 
effects at the cluster boundary may affect the results. In 
the PI model, the "cluster" is reduced to its minimum 
size, a single ion, and cluster-lattice interactions are de- 
scribed selfconsistently in the framework of the TESIm 
The cluster approximation can be rigourously formulated 
within the TES, as cluster-lattice orthogonality is a fun- 
damental requirement of that theory. Another feature of 
the PI model is that it does not invoke the LCAO approx- 
imation. The one-center character of the model leads to 
a large computational saving with respect to any multi- 
center cluster approach. The PI model just described has 
been formulated at the Hartree-Fock level. Now we turn 
to the introduction of Coulomb correlation. 



B. Introduction of correlation in the PI model. The 
Coulomb hole model 

The seminal idea of separability is to find electronic 
groups such that intergroup correlations play a minor 



roleSi Intraatomic correlation, on the other hand, con- 
tributes significantly to the band gap energy of ionic crys- 
tals and has to be included. To take into account in- 
traatomic correlation we have used the Cou]paab.-Hartree- 
Fock (CHF) model proposed by Clementi23H. In the 
CHF model, the coulomb-repulsion integrals are modi- 
fied by introducing a spherical hole around each electron 
in which the other electrons do not penetrate. The radius 
of this coulomb hole for a particular integral depends on 
the overlaps between the functions involved, aside from 
other factors. Two parameters can scale the hole, and 
their values were chosen to match the empirical correla- 
tion energies of He and Ne. Details <j*f-i|the CHF model 
can be consulted in dementi's paperst3H. 

From the point of view of the PI calculations, the 
CHF model is convenient for several reasons: (a) it is 
computationally simple; (b) it can be implemented eas- 
ily in the_B,oothaan-Bagus SCF (Self Consistent Field) 
formalisnfj; (c) it reproduces within a few parts per cent 
the empirical correlation energy for the ground state of 
free atomsoo; (d) it is size consistent; and (e) the CHF 
correlation energy depends on the radial density of the 
ions, thus incorporating correlation effects upon the for- 
mation of the solid. 

There are two ways to implement the CHF model. 
Without going into more detail, we have used the un- 
relaxed CHF approximation, in which the PI wavefunc- 
tions calculated at the HF level are kept fixed and the 
correlation energy is simply added to the HF energies. 
Since the coulomb repulsion integrals are reduced by the 
effect of the hole, the corrected energies of the ions are 
smaller than their corresponding original PI values. The 
good performance of the CHF model in accounting for in- 
traatoinic correlation can be appreciated in the reviews 
of refskao 



C. Basis set 

The localized atomic-like orbitals used to describe each 
ion in. the , crystal are expanded in a large basis set of 
STO'sEHtia because of their superior performance (their 
precision being near Numerical-Hartree-Fock). An ex- 
tension of the basis was required for the description of 
bromides. From an analysis of preliminary calculations 
we have reasons to think that the basis set used for Br~ 
may not be the most appropriate to describe this anion in 
the crystalline environment. Here we have decided to en- 
large the basis set for Br~ by adding one 5p-polarization 
function. The exponent of the 5p-STO was adjusted by 
minimizing the total energy of the corresponding crystal. 
We have found that upon inclusion of the 5p-STO, the 
4p-orbital of Br~ experiences a small contraction with 
respect to the 4p-orbital in the absence of the polariza- 
tion function. This very small effect is responsible for the 
improvement of the crystal energy. Within the TES, the 
portion of the total energy represented by the projection 
energy is sensitive to the quality of the basis set in the 
tail region. Further discussion on the influence of basis 
sets on the results is provided in section III. 
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III. CALCULATION OF ENERGY GAPS 

Following the work of Poole et a/Jil, the gap is rig- 
orously defined as the difference between the threshold 
energy E t and the electronic affinity \ °f the cristal: 

AE gap =E t - X - (13) 

The threshold energy is the energy needed to remove an 
electron from the top of the valence band (VB), whereas 
the electronic affinity is defined as the energy of an elec- 
tron at the bottom of the conduction band (CB), referred 
to the vacuum level, ft is important to be aware of this 
definition, because if the hole in the VB and the electron 
in the CB were allowed to interact (electron-hole pair), 
strictly we would be describing instead an excitonic level. 

The problem for a cluster model when trying to ob- 
tain band gap energies using eq.(13) is the calculation of 
the electronic affinity x, as it is necessary to deal with 
an electron in a [delocalized state at the bottom of the 
conduction bandit. The PI model is not an exception 
in this respect. Due to the one-center character of the 
model, we have to center the electron wave function on 
a lattice site. Besides, the strong ion-lattice orthogonal- 
ity required by the TES would force the electron to be 
localized on a given ion, leading to an incorrect represen- 
tation of a conduction band state. We can, however, give 
an approximate description of this state by relaxing the 
orthogonality requirements (see below). For ionic mate- 
rials the experimental values of \ are small (only a few 
tenths of eVE3o) compared to E t and our approximation 
gives values for x within the correct order of magnitude. 
Besides, typical errors in measured gaps are ~ 0.5 eV. 

We have calculated the band gap energy in two dif- 
ferent ways. The first one identifies the gap with the 
difference between the energy eigenvalues corresponding 
to the lowest unoccupied molecular orbital (LUMO) and 
the highest occupied molecular orbital (HOMO) obtained 
in the PI model at the HF level, that is, we approximate 
the threshold energy E t by the HOMO and the electronic 
affinity \ by the LUMO (with opposite signs): 

AE gap = e{LU MO) -e(HOMO). (14) 

This energy difference overestimates the_gap, as is the 
case of typical band theory calculationsElJ. The second 
way is a ASCF calculation. 

A. LUMO - HOMO difference. 

The orbital energy of the HOMO level is a quantity 
readily obtained in the PI model. For a given alkali halidc 
crystal it corresponds to the eigenvalue of the outermost 
occupied p-orbital of the anion. To obtain the LUMO we 
can simulate a neutral alkali atom, ^4° , as an impurity in 
the field created by the pure crystal (A : A + X~). Then, 
the LUMO is identified with the outermost occupied or- 
bital of the neutral alkali impurity. As stated above, 
when the calculation is performed within the strict frame- 
work of the PI model, the localized basis and the condi- 
tion of strong orthogonalization between the orbitals in 



neighboring sites would lead to an unphysical localization 
of this electron on the alkali site. A better description of 
the LUMO is achieved by freezing the crystal around the 
A impurity (that is the wave functions of the surround- 
ing ions) and removing the pieces Vx(i) and P s (i) out 
of the effective Hamiltonian for the impurity. It has to 
be stressed that our calculation describes the LUMO as 
an occupied orbital centered on an alkali atom. To allow 
for the derealization of this orbital over a substantial 
region of the crystal, the basis set .used for the impurity 
(taken from Bung©, and Barrientoscd for Li, Na, K, Rb, 
and from McLeanES for Cs) has been enlarged by adding 
some diffuse s-type STO functions, assuming that their 
exponents form a geometric progression^ For the ratio 
of the geometric progression we took the ratio between 
the exponents of the two outermost s-orbitals of the orig- 
inal basis set, and we enlarged successively the basis set 
with one, two and three STO's. The calculated LUMO 
eigenvalue converges fast as the basis set is enlarged. In 
fact, by adding a fourth diffuse function, the LUMO en- 
ergy changed by ~ 10 -3 eV only. With the basis chosen, 
the electron in the LUMO orbital of alkali halides spreads 
over a range of about 32 coordination cells from its cen- 
ter. By removing the pieces Vx(i) and P s (i) from the 
effective hamiltonian for the impurity, our calculation of 
the LUMO becomes rather similar to the method used in 
some embedded cluster models, in which the metal im- 
purity, with an extra electron added, is surrounded by a 
set of point charges chosen and optimized to represent 
the Madelung potentialc 3 ! The difference between the 
two procedures is that in our calculation we employ the 
complete charge distribution of the crystal and not ionic 
point charges. 

The calculation is similar for the oxide crystals. 
£homo is the eigenvalue of the highest occupied anionic 
orbital in the pure crystal, and (.lumo is obtained by 
simulating a metal (or silicon) cation with an electron 
added in its external s shell as an impurity in the pure 
crystal. Unlike the case of alkali halides, in the oxides the 
metal (or Si) cation remains charged (Be + , Mg + , Ca + , 
Sr+, Ba+, Al 2+ , Sift,) after addition of that electron. 
The Clementi-RoettiEU basis for the cations, enlarged by 
some diffuse functions, were used to simulate the der- 
ealization of the LUMO orbital. 

The magnitude of the LUMO eigenvalue is always a 
small quantity (the largest value is 0.4 eV for Si 3+ : 
SiO^)- Our results concerning the affinity lejffil agree 
with those obtained for MgO by Bagus et alSB in the 
sense that this level bears no simple connection with 
cationic s-orbitals. In all cases it was irrelevant to center 
the electron wave function on a cationic or an anionic 
site, the only important issue was to allow for the neces- 
sary derealization of the electron. The number of point- 
charge shells included in our calculation to converge the 
Madelung term is about 190, and the picture that arises 
is that the electron becomes as delocalized-as possible, 
the same conclusion reached by Bagus et atB. 

Results for the calculated band gaps of twenty al- 
kali halid e , . crvst als are compared in Table I with 
experimentu'c3c.a and with other calculations. In a 
similar way, Table II shows the results for seven ox- 
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ide crystals, also comparjed. with experimentc3ca and 
other theoretical resultsOCJ. .Typical errors in the ex- 
perimental gaps are 0.2-0.5 eVcll. Our calculations were 
performed using the experimental geometries and lat- 
tice constants of the perfect crystalsEj'Eil. The calcu- 
lated (LUMO — HOMO) gaps display the experimental 
trends: a decrease of the band gap by moving down the 
periodic table along the cation column when the anion is 
fixed (the experimental gaps for Br and I-crystals show 
some exceptions to this trend), and also a decrease by 
moving down the anion column when the cation is fixed. 
For the alkali halides, we obtain a general overestima- 
tion of the gap, with errors which increase with the size 
of the ions involved. The same trends are observed in a 
calculation by Bcrrondo and Rivas-Silvat3 using a model 
of a cluster embedded in a background of point-like ions. 
Berrondo's calculation was performed at the HF level 
and the gap, obtained as the difference of LUMO and 
HOMO orbital energies, is given in Table I. The error in 
the present approximate calculation is, however, about 1 
eV smaller than in Berrondo's calculation. On the other 
hand, for oxide crystals, the {LUMO — HOMO) calcu- 
lation yields excellent agreement with experiment. At 
this level of theory, correlation and orbital relaxation ef- 
fects are not yet taken into account, so this exceptional 
agreement has to be considered as resulting from a subtle 
cancellation of errors. 



B. ASCF calculation. 

The energy gap of an N-electron system is rigorously 
defined within the ASCF method as 

AE gap — (En—i — En) — (En — En+i), (15) 

where En is the ground state energyH. This is, evidently, 
equivalent to eqn. (13). For the alkali halide crystals we 
can write 

E t = Ejf-i — En — 
E cr ystai{X° : A+X-) - E crystal (A + X-), (16) 

X — ^n — En+i = 
E crystal {A+ X-) - E crystal (A° : A+X~), (17) 

associated to the two independent processes of removing 
an electron from a halogen anion and of placing an elec- 
tron in a state at the bottom of the conduction band, 
respectively. E crysta i{X : A + X~) represents the en- 
ergy of the crystal with a single neutral halogen impurity 
and E cryst ai(A Q : A + X~) that of a crystal with a "neu- 
tralized" alkali ion (although the neutralizing electron is 
in a fully delocalized state). Finally E crysta i(A + X~) is 
the energy of the perfect crystal. The geometry of these 
"doped" crystals is that of the perfect crystal, according 
to the Franck-Condon principle. 

The electronic orbitals of the ions surrounding the neu- 
tral impurities (in fact, of all ions in the crystal) were kept 
frozen in the same electronic states as in the pure crystal. 
This is not strictly necessary for the case X° : A + X~ but 



we have found from test calculations that the electronic 
relaxation of the cations surrounding the neutral halogen 
atom is indeed negligible. For the case of A : A + X~ , 
freezing of the electronic clouds of the surrounding shells 
of ions is convenient in view of the treatment of the im- 
purity (see below). With these assumptions the energies 
E t and \ reduce to a difference of effective energiesE-3 

Et = E eff (X° : A+X-) - E eff (X~ : A + X~) (18) 

X = E eff (A+ : A+X-) - E eff (A° : A+X~) (19) 

where the effective energies of J" and A + refer, evi- 
dently, to ions of the perfect crystal. Again, as in the dis- 
cussion of the HOMO-LUMO gap in section III. A above, 
we have enlarged the basis set to be able to describe 
an electronic state formally centered on an alkali site al- 
though delocalized over a region corresponding to many 
shells of neighbours. Applying the same argument to 
alkaline-earth oxides, the equivalent of equations (18) and 
(19) are 

E t = E eff (0- : A 2+ 2 -) - E eff (0 2 - : A 2+ 2 ~) (20) 

X = E eff (A 2+ : A 2+ 2 -) - E eff (A+ : A 2+ 2 -), (21) 

and related expressions can be easily written for the S1O2 
and AI2O3 crystals. 

The main contribution to AE gap comes from E t , and 
X only provides a small correction. The largest value of 
X, which is ~ 0.5 eV is found for Si02- This magnitude 
is consistent with the energy difference between the vac- 
uum level and the lowest conduction band state measured 
by photoemission spectroscopy^]. That energy difference 
is generally smaller than 1 eV, and in most cases much 
lower. 

The gaps obtained by this ASCF calculation are given 
in tables I and II. Two sets of results are included: the 
first one corresponds to calculations at the Hartree-Fock 
level, and the second includes correlation via the CHF 
model. The gaps obtained for alkali halides are smaller 
than those calculated by substracting LUMO and HOMO 
eigenvalues, a fact that improves the agreement with ex- 
periment, in particular for cases involving heavy halo- 
gens. The improvement arises from orbital relaxation 
in response to the removal of one electron from the an- 
ion. The relaxation of the orbitals of the halogen atom 
is included now in the calculation of E tl but not in the 
calculation of the HOMO level in section III. A. This ef- 
fect of a higher sensitivity of the anions (as compared to 
cations) to orbital relaxation-effects was also found in our 
previous works on clustersLru. For oxides, improvement 
with respect to the LUMO — HOMO difference of eqn. 
(14) is only found after including correlation. 

The electron density of the cations is very localized 
and practically does not change when the delocalized 
electron is centered on those cations. Furthermore, the 
delocalized electron does not interact with itself. This 
means that correlation gives a negligible contribution to 
X, so one can concentrate on the contribution of corre- 
lation to E t . The difference between ASCF(HF) and 
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ASCF(CHF) in Tables I and II comes from this source. 
We can appreciate in Table I that the correlation correc- 
tion improves the gap for most alkali halides, with KF, 
KC1, RbCl and CsCl as the only exceptions. Further- 
more, it is worth noting that this correction changes the 
magnitude of the gaps in the correct direction, namely, 
increasing the gap energy for fluorides, which have HF 
gaps below the experimental values, and decreasing the 
gap for the rest of the crystals, which show (except KC1, 
RbCl and CsCl) HF gaps larger than the experimental 
ones. In summary, the global trends remain the same as 
in the HF calculations, but the quantitative agreement 
improves with the inclusion of correlation effects. The 
same can be said of the oxide crystals, where inclusion of 
correlation reinstates the initial agreement found in the 
(LUMO — HOMO) calculation. These results lead us to 
conclude that correlation corrections are important for a 
quantitative comparison with experiment, although the 
general trends are reproduced at the J3F level. 

The band gaps obtained by Kunzp. for twelve alkali 
halide crystals and by Pandey et al23 for three oxides 
from band structure calculations are given for compari- 
son in Tables I and II. Energy band calculations at the 
HF level OApsestimate the band gap with errors of a few 
eV (see refHJ); notice.±hat this is a LUMO-HOMO-type 
of estimation. Kunzcirnij then introduced coulomb corre- 
lation by updating a method described by Pantelides et 
al£3. Correlation effects lowered the magnitude of the 
gap, bringing the results (included in Tables I and II) in 
good agreement with experiment. Comparing with the 
results obtained by Kunz, we can see that our CHF re- 
sults are slightly better, except for a few crystals (NaF, 
KF, RbBr), and that Kunz calculations are less success- 
ful in reproducing the general trends. We have just men- 
tioned that HF band structure calculations overestimate 
the gap by a sizable amount. On the other hand our 
ASCF calculations at the HF level do not seem to be 
much affected by this problem. We ascribe this to elec- 
tronic relaxation effects. 

In Table J_ we also include the results of Berrondo and 
Rivas-Silvao, obtained as the difference between the to- 
tal energies of the lowest electronic excited state and 
the ground state of a cluster simulating a piece of the 
ionic crystal, embedded in a system of point charges. 
These calculations were performed at the HF level and, 
evidently, include electronic relaxation, so in Table I 
the corresponding results have been classified under the 
ASCF category. Berrondo's calculations include in addi- 
tion relativistic effects. A comparison with the present 
ASCF(HF) results shows rather similar gaps. In the 
same paper Berrondo and Rivas-Silva have obtained the 
gap by identifying this one with the excitation energy 
obtained in a Configuration Interaction (CI) calculation, 
again within a cluster model, in this case a bare cluster 
without surrounding point charges. As a consequence of 
the introduction of correlation, the CI calculations give 
an improvement of the calculated gaps for chlorides, bro- 
mides and iodides. In the case of fluorides, where some 
of the ASCF(HF) gaps obtained by Berrondo are below 
the experimental values, the CI calculation leads to an 
increase of the gaps. So turning from ASCF to CI seems 



to work in the right direction, although for fluorides this 
results in a worsening of the magnitude of the gap. The 
errors in the CI gaps are not larger than I eV. This is also 
the case for our ASCF calculations including correlation, 
except for NaF, KC1 and RbCl, for which the error is 1.1 
eV, and for CsCl and CsBr, where the errors are 1.6 eV 
and 1.9 eV respectively. The computer time required by 
our calculations is, however, much smaller. 

Included in Table II are the theoretical results obtained 
by Lorda et uo for alkaline-earth oxides. They used a 
model of a linear triatomic M-O-M cluster (M= metal, 
0= oxygen) embedded in a Madelung field representing 
the crystal. The ab initio wavefunction was written as a 
linear combination of Slater determinants, each one cor- 
responding to a resonating valence bond structure. In 
that way a valence bond wave function is obtained which 
is explicitely correlated from the outset, in the sense that 
it is not an eigenfunction of a hamiltonian which can be 
written as a sum of monpelectronic operators (for details, 
consult chapter 7 of refEl). In order to include external 
correlation effects they used standard second-order per- 
turbation techniques. This external correlation is defined 
as the energy contribution of those Slater determinants 
not explicitely included in the valence bond model space. 
The gap was then calculated as an energy difference be- 
tween appropriate resonating valence bond components. 
Our results are also in better agreement with experiment 
in this case. Actually, the first excitation energies calcu- 
lated by Lorda et al. should not be identified with the 
band gap energies, but with localized excitations. This 
point was emphasized when we gave the definition of the 
band gap energy above. ■— ■ 

Density functional theory (DFT)E3 has also been used 
to calculate band gaps in ionic crystals. The DFT one- 
particle eigenvalues have no formal justification as quasi- 
particle energies although in practice these eigenvalues 
have been used to discuss the spectra of solids. The DFT 
band gaps (usually calculated employing the local density 
approximation for exchange and correlation effects) give 
values typically 30-50 % less than the gaps observed in 
the optical spectra. A rigorous formulation requires the 
calculation of the quasi-particle band-structure. This has 
been done by Hybertsen and LouieO by usuk the so- 
called GW approximation to the self-energjC3. These 
authors have obtained a gap of 9.1 eV for LiCl. 

IV. CONCLUSIONS. 

A preliminary requirement for the treatment of lumi- 
nescent impurities in ionic crystals is an accurate descrip- 
tion of the band gap. Consequently, we have carried out 
a series of approximate calculations of the band gap for 
twenty alkali halide and seven oxide crystals. Two differ- 
ent prescriptions have been used: the first one is a differ- 
ence between the energies of the LUMO and HOMO or- 
bitals of the crystal, and the second is a ASCF calculation 
based on subtracting total energies. For this purpose the 
Perturbed Ion (PI) model has been used, supplemented 
with some additional assumptions to deal with the state 
corresponding to one electron at the bottom of the con- 
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duction band. Correlation effects have been taken into 
account by a model introduced by Clementi. 

Overall we have obtained a good description of the en- 
ergy gaps with a very modest computational effort. The 
gaps show a reasonable agreement with experiment, and 
also compare well (in many cases favorably) with previ- 
ous theoretical works. The results support the validity of 
the PI model and the additional assumptions employed 
to simulate the physical situation. Further improvement 
of some particular points (as the modelling of the delo- 
calized electron at the bottom of the conduction band, 
the treatment of electronic correlation, and the introduc- 
tion of relativistic effects for heavy atoms) are expected 
to improve further the results. 



Acknowledgments: This work has been supported 
by DGICYT(Grant PB95-0720-C02-01). A. Aguado is 
greatful to Univcrsidad de Valladolid, Caja Salamanca 
y Soria, and Junta de Castilla y Leon for Postgraduate 
fellowships. A. Ayuela thanks the hospitality and sup- 
port from Brigham Young University during a visit at 
the starting of this work. 



1 G. BLASSE and B. C. GRABMAIER, Luminescent mate- 
rials, Springer Verlag, Berlin (1995). 

2 M. BERRONDO and J. F. RIVAS-SILVA, Int. J. Quantum 
Chem. Quantum Chem. Symp. 29, 253 (1995). 

3 S. HUZINAGA and A. P. CANTU, J. Chem. Phys. 55, 
5543 (1971). 

4 S. HUZINAGA, D. MCWILLIAMS, and A. A. CANTU, 
Adv. Quantum Chem. 183 (1973). 

5 V. LUANA. Thesis Disertation, Universidad de Oviedo 
(1987). 

6 V. LUANA and L. PUEYO, Phys. Rev. B 41, 3800 (1990). 

7 A. AYUELA, J. M. LOPEZ, J. A. ALONSO, and V. 
LUANA, Z. Phys. D 26, S213 (1993); {ibid.) Physica B 
212, 329 (1995); {ibid.) Canad. J. Phys. (in press). 

8 A. AGUADO, A. AYUELA, J. M. LOPEZ, and J. A. 
ALONSO, J. Phys. Chem. 101B, 5944 (1997); Phys. Rev. 
B, 56, 15353 (1997); E. DE LA PUENTE, A. AGUADO, 
A. AYUELA, and J. M. LOPEZ, Phys. Rev. B 56, 7607 
(1997). 

9 T. D. CLARK and K. L. KLIEWER, Phys. Lett. 27 A, 167 
(1968). 

10 A. B. KUNZ, Phys. Rev. 175, 1147 (1968). 

11 C. Y. FONG and M. L. COHEN, Phys. Rev. Lett. 21, 22 

(1968) . 

12 A. B. KUNZ, Phys. Stat. Sol. (b) 29, 115 (1968). 

13 W. B. FOWLER and A. B. KUNZ, Phys. Rev. 186, 956 

(1969) . 

14 L. J. PAGE and E. H. HYGH, Phys. Rev. B 1, 3472 (1970). 

15 A. B. KUNZ, J. Phys. C: Solid St. Phys. 4, 1542 (1970). 

16 H. OVERHOF, Phys. Stat. Sol. (b) 43, 575 (1971). 

17 N. O. LIPARI and A. B. KUNZ, Phys. Rev. B3, 491 (1971). 

18 F. PERROT, Phys. Stat. Sol. (b) 52, 163 (1972). 

19 R. T. POOLE, J. LIESEGANG, R. C. G. LECKEY, and 
J. G. JENKINS, Phys. Rev. B 11, 5190 (1975). 



20 C. JOUANIN, J. P. ALBERT, and C. GOUT, II Nuovo 
Cimento 28, 483 (1975). 

21 A. B. KUNZ, Phys. Rev. B 26, 2056 (1982). 

22 R. A. HEATON, J. C. HARRISON, and C. C. LIN, Phys. 
Rev. B 28, 5992 (1983). 

23 R. PANDEY, J. E. JAFFE, and A. B. KUNZ, Phys. Rev. 
£43, 9228 (1991). 

24 P. G. LYKOS and R. G. PARR, J. Chem. Phys. 24, 1166 
(1956). 

25 R. G. PARR, R. O. ELLISON, and P. G. LYKOS, J. Chem. 
Phys. 24, 1106 (1956). 

26 S. HUZINAGA, L. SEIJO, Z. BARANDIARAN, and M. 
KLOBUKOWSKI, J. Chem. Phys. 86, 2132 (1987). 

27 H. J. SILVERSTONE and R. K. MOATS, Phys. Rev. A 
16, 1731 (1977). 

28 A. MARTIN PENDAS and E. FRANCISCO, Phys. Rev. 
A 43, 3384 (1991). 

29 Z. BARANDIARAN and L. PUEYO, J. Chem. Phys. 79, 
1926 (1983). 

30 A. B. KUNZ and J. M. VAIL, Phys. Rev. B38, 1058 (1988). 

31 R. MCWEENY, Methods of molecular quantum mechanics, 
Academic Press. London (1994). 

32 E. CLEMENTI, IBM J. Res. Develop., 9, 2 (1965). 

33 S. J. CHAKRAVORTY and E. CLEMENTI, Phys. Rev. A, 
39, 2290 (1989). 

34 C. C. J. ROOTHAAN and P. S. BAGUS, Methods Comput. 
Phys., 2, 47 (1963). 

35 L. PUEYO, V. LUANA, M. FLOREZ, and E. FRAN- 
CISCO in Structure, Interactions and Reactivity, edited by 
S. FRAGA (Elsevier, Amsterdam, 1992), vol. B, p. 504. 

36 V. LUANA, M. FLOREZ, E. FRANCISCO, A. MARTIN 
PENDAS, J. M. RECIO, M. BERMEJO, and L. PUEYO, 
in Cluster Models for Surface and Bulk Phenomena, edited 
by G. PACCHIONI, P S. BAGUS and F. PARMIGIANI 
(Plenum, New York, 1992), p. 605. 

37 E. CLEMENTI and C. ROETTI, At. Data Nucl. Data Ta- 
bles, 14, 3 (1974); 14, 177 (1974). 

38 A. D. MCLEAN and R. S. MCLEAN, At. Data Nucl. Data 
Tables, 26, 197 (1981). 

39 P. S. BAGUS, F. ILLAS, and C. SOUSA, J. Chem. Phys. 
100, 2943 (1994). 

40 R. T. POOLE, J. G. JENKINS, R. C. G. LECKEY, and 
J. LIESEGANG Phys. Rev. B 11, 5179 (1975). 

41 C. F. BUNGE and J. A. BARRIENTOS, At. Data Nucl. 
Data Tables 53, 113 (1993). 

42 R. C. RAFFENETTI and K. RUEDENBERG. Even- 
Tempered Representations of Atomic Self-Consistent-Field 
Wave Functions (Ames Laboratory Report, 1973). 

43 J. CASANOVAS, C. SOUSA, J. RUBIO, and F. ILLAS, J. 
Comput. Chem. 14, 680 (1993). 

44 W. HAYES and A. M. STONEHAM, Defects and Defect 
Processes in Non-Metallic Solids, Wiley, New York (1985). 

45 Handbook of Laser Science and Technology, edited by M. 
J. WEBER (CRC, Cleveland, 1986) Vol. III. 

46 W. H. STREHLOW and E. L. COOK, J. Phys. Chem., 
Ref. Data 2, 163 (1973). 

47 E. T. ARAKAWA and M. W. WILLIAMS, J. Phys. Chem. 
Solids 29, 735 (1968). 

48 F. J. HIMPSEL, Surf. Sci. 168, 764 (1986); F. J. GRUN- 
THONER and P. J. GRUNTHONER, Mater. Sci. Rep. 1, 
65 (1986). 

49 A. LORDA, F. ILLAS, J. RUBIO, and J. B. TORRANCE, 
Phys. Rev. B 47, 6207 (1993). 

50 R. W. G. WYCKOFF, Crystal Structures, Volume I, New 



7 



York: Wiley-Interscience (1963). 

C. KITTEL, Introduction to Solid State Physics, Wiley, 
New York (1976). 

M. BERRONDO and J. F. RIVAS-SILVA. Int. J. Quantum 
Chem., 57, 1115 (1996). 

R. W. GODBY, M. SCHLUTER, and L. J. SHAM, Phys. 
Rev. B 37, 10159 (1988); L. J. SHAM and M. SCHLUTER, 
Phys. Rev. B 32, 3883 (1985). 

A. MARTIN PENDAS, E. FRANCISCO, V. LUANA, and 
L.PUEYO, J. Phys. Chem. 96, 2301 (1992). 
S. T. PANTELIDES, D. J. MICKISH, and A. B. KUNZ, 
Phys. Rev. B 10, 2602 (1974). 

Theory of the Inhomogeneus Electron Gas. Editors S. 
LUNDQVIST and N. H. MARCH, Plenum, New York 
(1983) 

M. S. HYBERTSEN and S. G. LOUIE, Phys. Rev. B 34, 
5390 (1986). 

L. HEDIN, Phys. Rev. 139, A796 (1965); L. HEDIN and 
S. LUNDQVIST, Solid State Phys. 23, 1 (1969). 



Captions of tables 

Table I Band gaps of alkali halide crystals calculated 
by diffcrcnt||theOji'etical methods, compared to experimen- 
tal valuesESSti Ae ffF (LUMO - HOMO) is a difference 
of one-particle orbital energy eigenvalues. ASCF refers 
to a difference between the energies of the crystal in ap- 
propriate electronic states (see text). KunzEil performed 
band structure calculations for the solid and CI repm- 
sents a Configuration-Interaction cluster calculations. 
All energies are given in eV. 

Table II Band gaps of seven oxides calculated from 
several theoretical methods, compared to experimental 
values. Achf is a (LUMO — HOMO) difference of one- 
particle orbital energy eigenvalues. ASCF refers to a dif- 
ference between the energies of the crystal in appropriate 
electronic states (see text). Also included are the ab ini- 
tio band structure calculations of Pandey et aa and the 
ab-initio valence bond cluster model results of Lorda et 
aa All energies are given in eV. 
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